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Abstract 

We propose an approach to compute the boundary crossing probabilities for a class of dif- 
t-H ■ fusion processes which can be expressed as piecewise monotone (not necessarily one-to-one) 

functionals of a standard Brownian motion. This class includes many interesting processes 
in real applications, e.g., Ornstcin-Uhlenbeck, growth processes and geometric Brownian mo- 
tion with time dependent drift. This method applies to both one-sided and two-sided general 
' nonlinear boundaries, which may be discontinuous. Using this approach explicit formulas for 

7-H \ boundary crossing probabilities for certain nonlinear boundaries are obtained, which are useful 

^sO • in evaluation and comparison of various computational algorithms. Moreover, numerical com- 

putation can be easily done by Monte Carlo integration and the approximation errors for general 
boundaries are automatically calculated. Some numerical examples are presented. 
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1 Introduction 



Let X = {X t ,t > 0} be a diffusion process denned on a probability space (Q,A,P) and has state 
space either the real space M or a subinterval of it. In this paper, we are concerned with the 
following boundary crossing probability (BCP) 

P x (a,b,T) = P(a(t) <X t < 6(t),Vt € [0,T]) , (1) 

where T > is fixed, boundaries a(t) and b(t) are real functions satisfying a(t) < b(t) for all 
< t < T and a(0) < x < 6(0). 

The problem of boundary crossing probabilities, or the first passage time (FPT) distributions, 
has drawn tremendous amount of attention in many scientific disciplines. To mention a few, for 
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example, it arises in biology (Ricciardi et al, 1999), economics (Kramer, Ploberger and Alt 1988), 
engineering reliability (Ebrahimi 2005), epidemiology (Martin-L6f 1998, Tuckwell and Wan 2000 
and Startsev 2001), quantitative finance (Garrido 1989, Roberts and Shortland 1997, Lin 1998, 
Novikov et al, 2003, Borovkov and Novikov 2005), computational genetics (Dupuis and Siegmund 
2000), seismology (Michael 2005), and statistics (Doob 1949, Anderson 1960, Durbin 1971, Sen 
1981, Siegmund 1986, Bischoff et al 2003, and Zeileis 2004). 

Despite its importance and wide applications, explicit analytic solutions to boundary crossing 
problems do not exist, except for very few instances. Traditionally, the mainstream of the research of 
nonlinear boundary problems is based on Kolmogorov partial differential equations for the transition 
density function, and focuses on approximate solutions of certain integral or differential equations 
for the first-passage time densities. For example, to obtain the approximate solutions for Brownian 
motion crossing continuously differentiable boundaries, the tangent approximation and other image 
methods have been used by Strassen (1967), Daniels (1969, 1996), Ferebee (1982) and Lerche 
(1986), whereas a series expansion method has been used by Durbin (1971, 1992), Ferebee (1983), 
Ricciardi et al. (1984), Giorno et al. (1989), and Sacerdote and Tomassetti (1996) to deal with more 
general diffusion processes. Recently, Monte Carlo path-simulation methods have been proposed to 
numerically compute the FPT density for a general diffusion process crossing a one-sided constant 
boundary (Kloeden and Platen 1992, Giraudo and Sacerdote 1999, and Giraudo, Sacerdote and 
Zucca 2001). 

In principle, all these FPT methods are designed for the continuously differentiable boundaries 
only, and they mostly deal with one boundary problems. However, two-sided and discontinuous 
boundary problems arise in many real applications, e.g., in problems of pricing barrier options in 
quantitative finance. Furthermore, the actual numerical computations in FPT methods are either 
intractable or give approximate solutions for which the approximation accuracies are difficult to 
assess (Daniels 1996, Sacerdote and Tomassetti 1996). The more recent path-simulation methods 
require heavy computation and the algorithms are typically very complicated (Giraudo, Sacerdote 
and Zucca 2001). 

In contrast to the traditional FPT methods, Wang and Potzelberger (1997) proposed an alter- 
native approach, which focuses on the boundary crossing probabilities (BCP) directly. Using this 
approach, Wang and Potzelberger (1997) derived an explicit formula for the BCP for Brownian 
motion crossing a piecewise linear boundary, and then used this formula to obtain approximations 
of the BCP for more general nonlinear boundaries. The numerical computation can be easily done 
by Monte Carlo integration, and the approximation errors are automatically computed. This ap- 
proach also allows one to control the approximation error before the actual computation begins. 
This method has been extended to two-boundary problems by Novikov, Frishling and Kordzakhia 
(1999), who use recurrent numerical integration with Gaussian quadrature, and by Potzelberger 
and Wang (2001), who use Monte Carlo integration. Moreover, it has been applied in finance 
(Novikov, Frishling and Kordzakhia 2003, Borovkov and Novikov 2005), change-point problems 
(Zeileis 2004), and engineering reliability (Ebrahimi 2005). 

In this paper, we extend the method of Wang and Potzelberger (1997) and Potzelberger and 
Wang (2001) further to a class of diffusion processes, which can be expressed as piecewise monotone 
(not necessarily one-to-one) functionals of a Brownian motion. This class contains many interesting 
processes arising in real applications, e.g., Ornstein-Uhlenbeck processes, growth processes and 
the geometric Brownian motion with time-dependent drift. Moreover, this approach allows us to 
derive explicit formulas for BCP for certain nonlinear boundaries, which are useful in evaluation 
and comparison of various computational methods and algorithms for boundary crossing problems. 

From methodology point of view, the framework used in this paper is different from that in the 
previous works in the literature. Whereas the previous works are based on the Kolmogorov partial 
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differential equations for the transition density functions, our approach is based on the stochastic 
differential equation (SDE) for the processes. The SDE approach becomes more and more popular 
in modern times, especially in mathematical finance since the publication of the stimulating works 
of Black and Scholes (1973) and Merton (1973). See, e.g., Garrido (1989), Roberts and Shortland 
(1997), Tuckwell and Wan (2000), Kou and Wang (2003), Novikov, Frishling and Kordzakhia (2003), 
and Ebrahimi (2005). 

The paper is organized as follows. Section 2 extends the BCP formulas of Wang and Potzelberger 
(1997) and Potzelberger and Wang (2001) for Brownian motion to piecewise continuous boundaries. 
Section 3 derives the BCP formulas for more general diffusion processes. Sections 4, 5 and 6 treat 
the Ornstein-Uhlenbeck processes, growth processes and geometric Brownian motion respectively. 
Finally, section 7 gives numerical examples and section 8 contains conclusions and discussion. 



2 Brownian motion and piecewise continuous boundaries 



Let W = {Wt,t > 0} be a standard Brownian motion (Wiener process) with EWt = and 
EW t W s = min(t,s). We consider the BCP P w (a,b) = P (a(t) <W t < b(t),\ft £ [0,T]). Wang 
and Potzelberger (1997), Novikov, Frishling and Kordzakhia (1999), and Potzelberger and Wang 
(2001) have derived explicit formulas for P\v(a,b,T) for continuous boundaries a and b. It is 
straightforward to generalize their results to piecewise continuous boundaries. Throughout the 
paper we denote by (ti)f =1 , < t\ < ■ ■ ■ < t n -\ < t n = T, a partition of interval [0, T] of size 
n > 1. Further, let to = 0, Aij = U — U-i, on = a(ti), 0i = b(U) and 5i = — on. Then completely 
analog to Wang and Potzelberger (1997) and Potzelberger and Wang (2001), we can establish the 
following results. 

Theorem 1 Given any partition (U)2 =1 of interval [0, T], suppose the boundaries a andb are linear 
functions in each subinterval (U-i,ti) and are such that, at any ti, lim^ ti _ a{t) > lim^ ti+ a{t) 
and \\m. t ^ti- b(t) < lim t ^. ti+ b{t) hold. Then 

P w (a,b,T) = Eg(W tl ,W t2 ,...,W t 

where function g{x), x = (xi,X2, ...,x n )' , is defined as follows. 

(1) For one-sided BCP P w (-oo, b,T), 

A r r 2 

g(x) = [[l { _ ooA) (x i )n-ex.p - — 
i=i 1 L ^ r 

where l(-) is the indicator function. 

(2) For two-sided BCP P w (a,b,T), 



■(Pi-i - Xi-i)([3i - Xi) 



(2) 



(3) 



g(x) = n i (a i ,/3 i )( a; <) 



i=l 



1 - ^2,hij{xi-i, 



J, ; 



(4) 



where 



hij (Xj_i , Xj) 



exp 



"At- + ^ ai ~ l ~ Xi ~ 1 ^ + ( ai ~ Xi 



exp | [jSi-idi + 6i-i(ai - Xi) - 5i(ai-i - Xj_i)]| 

exp {~~AT- ~ ~ Xi -^ ti Si ~ (^» ~ Xi )l} 
exp j-^: [jSi-iSi - 5i-i(Pi - Xi) + 6i(Pi-i - . 
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Note that function hij(xi-i,Xi) in ^) consists of exponential functions that decrease very rapidly 
with the increase of j. In all of our numerical examples, fairly accurate approximations are achieved 
by using six terms only. Moreover, a simplified version of is given in Potzelberger and Wang 
(2001). 

The results of Theorem [1] can be used to approximate the BCP Pw(a,b) for general nonlinear 
boundaries a and b, provided they can be sufficiently well approximated by some piecewise linear 
functions. More precisely, if the sequence of piecewise linear functions a n (t) — > a(t) and b n (t) — * b(t) 
uniformly on [0,T], then it follows from the continuity property of probability measure that 

lim P w (a n , b n ,T) = P w (a, b, T). (5) 

n — >oo 

The accuracy of approximation ([SJ depends on the partition size n. In general, larger n will give 
more accurate approximation. For twice continuously differentiable boundaries, Novikov, Frishling 
and Kordzakhia (1999) showed that the approximation errors converge to zero at rate \/log n/n 3 , 
if equally-spaced piecewise linear approximation is used. Later, Potzelberger and Wang (2001) 
obtained an approximation rate of 0(l/n 2 ) by using an "optimal" partition of [0, T]. More recently, 
Borovkov and Novikov (2005) obtained the same rate by using equally-spaced partition. 

In practice, one simple and straightforward way to assess the approximation accuracy is to 
calculate the lower and upper bounds of the BCP using two piecewise linear functions approaching 
a and b from each side respectively. This method will be demonstrated through numerical examples 
in section 7. See Wang and Potzelberger (1997) and Potzelberger and Wang (2001) for more details 
and examples. 

3 General diffusion processes 

In this and the subsequent sections, we generalize the results of the previous section to more general 
diffusion processes satisfying the following stochastic differential equation (SDE) 

dX t = fi(t, X t )dt + a(t, X t )dW t , X = s , (6) 

where the drift fj,(t, x) : [0, oo) x M — > M and diffusion coefficient a(t, x) : [0, oo) x M — > M + are real, 
deterministic functions and {Wt,t > 0} is the standard Brownian motion (BM) process. In order 
for the SDE © to admit unique solution, throughout this paper we assume generally that fi(t, x) 
and a(t,x) are measurable functions and satisfy the following Lipschitz and growth conditions: for 
any t > 0, 

\n(t,x) - fi(t,y)\ + \a(t,x) - a(t,y)\ < L\x-y\ ,Vx,y € M, (7) 

and for any x 6 JR, 

\n(t,x)\ + |<7(t,x)| < K{\ + 1x1), Vt > 0, (8) 

where L > and K > are constants. However, since the solution problem is not a major concern 
of this paper, in the subsequent derivations we will not deal with conditions ((7J) and flSJ) explicitly. 
Rather we refer the reader to, e.g., Karatzas and Shreve (1991, p289) or Lamberton and Lapeyre 
(1996, p50). See also Yamada and Watanabe (1971). In addition, throughout the paper we make 
the following assumption, which is necessary for the application of Ito's formula. 

Assumption 1 Let {Xt,t £ [0,T]} denote the solution of the SDE {2P with state space an open or 
closed finite or infinite interval S. Assume that a(t,x) > for all t 6 [0,T] and x £ S. Further, 
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assume that there exists a function f(t,x) G C 1,2 ([0,T] x S) (once continuously differentiate with 
respect to t and twice with respect to x), such that 1. for x G S 

df(t,x) = 1 
dx a(t, x) 

and 2. Ito's Lemma may be applied to Y t = f(t,X t ), i.e. 

Now we consider the BCP ([T]) for any diffusion process X = {X t ,t > 0} satisfying §§§ and 
Assumption [TJ The basic idea is that Px (a, b, T) can be computed through the BCP for the 
standard BM, as long as X can be expressed as a "nice" transformation of it. More precisely, given 
any boundaries a(t), b(t) and T > 0, let us denote 

A(a,b,T) = {x(t) G C([0,oo)) : a(t) < x{t) < b{t),Vt G [0,T]}. 

Then it is easy to see that, if there exists a standard BM W = {Wt,t > 0} and a measurable 
functional /, such that X = fiW), then 

P(X G A(a, b, T)) = P{W G f' l A{a, b, T)). (10) 

Consequently, the BCP for X can be calculated through the BCP for W, as long as there exist 
boundaries c(t), d(t) and S > 0, such that f~ 1 A(a,b,T) = A(c,d,S). One of such instances is 
described in the following theorem. 

Theorem 2 If there exists a function f(t,x) G C 1,2 ([0,oo) x M), such that Y t := f(t,X t ) satisfies 
dY t = d t dWt, where at G C([0, oo)) is a real, deterministic function satisfying &t ^ 0, Vi G [0, oo), 
then there exists a standard BM {W s , s > 0}, such that for any boundaries a(t),b(t), 

1. P x (a,b,T) = P w (c,d,S), ifa t >0, Vi G [0, oo); 

2. P x (a,b,T) = P w (d,c,S), ifa t <0, Vi G [0, oo); 
where 

c(s) = f(t(s),a(t(s))) - /(0, so), < s < S, (11) 
d(s) = f(t(s), b(t(s))) - /(0, x ), < s < S, (12) 
t(s) is the inverse function of s(t) = f^a^du and S = s(T). 

Proof. First, from dYt = dtdWt it follows that {Yt,t > 0} is a continuous Gaussian process with 
variance function ^ 

V(Y t ) = [ a 2 u du = s(t),t>0. 
Jo 

Because a\ > 0, s(t) is strictly increasing for t > and s(0) = 0. For every s > 0, define 
W s := Y t r s \ — Yq, where t(s) is the inverse function of s(t). Then since Wq = and 

V(W s ) = V(Y t{s) ) = s(t(s)) = s, 

{W s , s > 0} is a standard BM. Furthermore, for every t > 0, Yt = Yq + W s m. 
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On the other hand, Assumption Q] and Ito's formula ([9]) imply that 

/ sdf ~ 

(T{t,X)— = <Jf 

Note that for every t > 0, 

Of = 5-t 
dx a(t, x) 

is continuous in x. It follows that, for any xq and x in the state space of {X t }, 

f(t, x) = a t [ X -j^du + h(t), (13) 
J xo a[t,u) 

where h(t) is an arbitrary function which does not depend on x. Further, since at ^ and 
a(t,x) > 0, for every t > 0, f(t,x) is a monotone function of x. For every t > 0, let g(t,x) be the 
inverse function of f(t, x). Then X t = g(t, Y t ). 

Now we consider case 1, where at > 0, Vt G [0, oo). In this case function f(t,x) is strictly 
increasing in x. It follows that 

P(a(t)<X t <b(t),Vte[0,T\) = p(a(t)<g(t,Yo + W a{t ))<Kt),Vte[Q,T\) 

= P (f(t, a(t)) <Y + W s{t) < f(t, b(t)),Vt e [0, T]j 
= P (c(s) <W S < d(s), Vs G [0, 5]) , 

where c(s) and d(s) are given by (fTTj) and (jT2j) respectively, and S = s(T). 

Similarly, case 2 follows from the fact that function f(t,x) is strictly decreasing in x. □ 



Remark 3 The above derivation is readily generalized to more complicated case where at has finite 
number of zeros. Without loss of generality, suppose at has k zeros < t\ < ti < ■ ■ ■ < t\~ < T. 
Then at changes signs over intervals i = 1,2,..., k + 1, where to = and tk+i = T. 

Denote by T + the union of all intervals in which at > 0, and by T_ the union of all intervals in 
which a t < 0. Then by (|13|) function f(t,x) is strictly increasing for t £ T + and decreasing for 
t G TL. Furthermore, let S + = s(T + ), S- = s(T_), where s(t) = J a\du. Then for any boundaries 
a(t),b(t), we have Px(a,b,T) = Pyy(c, d, S), where 



c(s) 
d(s) 



f{t(s),a(t(s))) - /(0, xq), when a G 5 4 

f(t{s),b(t{8)))-f(p,x ), when seS, 

f(t(s),b(t(s)))-f(0,x ), when S G5 4 

f(t(s),a(t(s))) - /(0, x ), when s G 5_ 



and S = s(T). 



There remains the question that for which processes {Xt} do functions f(t, x) and at in Theorem 
[2] exist. A class of such diffusion processes is characterized by the next Theorem. 

Theorem 4 Let {Xt, t > 0} satisfy equation (0) with drift fj,(t, x) and diffusion coefficient a(t, x) G 
C 1,2 ([0, oo) x M). Then there exist functions f(t,x) and at satisfying conditions of Theorem^ if 
and only if lit = ii{t, x) and at = a(t, x) satisfy the following partial differential equation 



d_ 

dx 



1 dot ^ d_ flda t _ fH 
at dt dx \2 dx at 



0. (14) 
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Proof. Suppose Y t = f(t, X t ) and dY t = &tdWt, where f(t, x) and at satisfy conditions in Theorem 
[2j Then Ito's formula Q implies 

df df a 2 d 2 f n 

and 

(16) 



From (|16p . we have 

d 2 f _ a t da t 
dx 2 a 2 dx ' 

Substituting the above equation into (|15p. we obtain 

df_ _ - flda t _ fH 
dt \2 dx at 

It follows that 

'«•*>- jf*-Gtr-£)* + *M- (17> 

where <p(x) G C 2 (1R). Differentiating (|17p with respect to x and using equation (|16p . we have 

— a (-— - f^i] du + = — (18) 

dx Jo u \2 dx a u ) dx at 

Further, differentiating (fT8|) with respect to t yields 

d (lda t m\ d a t 

which is equivalent to 



at dx V 2 dx at ) dtat 



1 dat d flda t fXt\ 1 dat 
at dt 0t dx V 2 dx at / a t dt 1 

Since the right-hand side of (|19p does not depend on x, differentiating both sides of this equation 
with respect to x yields dHJ). 

Conversely, if ([T3| holds, then equation (fTUj) defines a function <7 t G C 1 ([0, oo)). Further, 
equation (fT8|) defines a function 0(x) G C 2 (M), whereas (fTT|) defines function f(t,x). From (fT7|) 
and (I14p . it is easy to see that f(t,x) satisfies condition in Theorem [2j □ 

Condition (|14p is practical in real applications. Once a process is given, one can easily check 
whether condition (|14p holds or not. On the other hand, equation (|14p can also be used to char- 
acterize certain subclasses of diffusion processes, which may be of interest of researchers working 
on particular problems. For example, for any given type of diffusion coefficient at, one can find all 
possible functions /it, such that (I14p is satisfied. Some examples are given in the following Corollary. 



Corollary 5 For each of the following types of diffusion coefficient at, the corresponding general 
solutions of equation f7^| ) for the drift fit are given as follows, where a{t), [3(t) are arbitrary real 
functions and a(t) : [0, oo) — ► M + . 

1. (L-class). For a(t, x) = a(t), fi(t, x) = a(t)x + (3(t), x G M. 

2. (G-class). For a(t, x) = a(t)x, fi(t, x) = a(t)x + j3(t)x log x, x G 
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The proof of these results is straightforward and is therefore omitted. It is easy to see that L-class 
corresponds to general Ornstein-Uhlenbeck processes with time-dependent coefficients, whereas G- 
class corresponds to growth processes which are widely used in population genetics. It is easy to 
see that the geometric BM process also belongs to G-class. 

There remains a practical question as how to find transformations f(t,x) and at for a given 
process. Obviously this depends on the functional forms of rtt and at and can only be worked out 
explicitly on case by case bases. Generally speaking, equation (fT9l) defines at, up to a multiplicative 
constant, and (|13j) defines f(t,x) up a term h{t) not depending on x, which is in turn determined 
by (fl5l) . In the subsequent sections, we demonstrate this for some widely used special processes 
and use the theoretical results of this section to derive more detailed results for these processes. 

Remark 6 The problem of obtaining the transition probability of a diffusion process through 
one-to-one (strictly monotone) transformation of Brownian motion was first posed by Kolmogorov 
(1931). A class of such diffusion processes was characterized by Cherkasov (1957) and Ricciardi 
(1976), whereas the most general class was derived by Bluman (1980). Using this transformation 
method, Ricciardi, Sacerdote and Sato (1984) obtained the FPT densities for a class of diffusion 
processes crossing one-sided, continuously differentiable boundaries. Since the main concern of this 
paper is the boundary crossing probabilities, not the transformation itself, we do not intend to be 
most general. The class of processes characterized by Theorem [J] is large enough to contain many 
interesting processes arising in practice. Furthermore, the condition (|14p is practical and easy to 
work with in real applications. 

4 Ornstein-Uhlenbeck processes 

Ornstein-Uhlenbeck (O-U) processes are a class of important diffusion processes with wide appli- 
cations. We start with the O-U process which is defined in state space M and satisfies 



where K,a G iR + and a G M are constants. In mathematical finance, (I20p is known as Vasicek 
model for the short-term interest rate process (Vasicek 1977). 

In this case, special forms of transformation f(t, x) and boundaries (llip and (jl2l) can be derived. 
For every t > 0, if we define 



dX t = n(a — X t )dt + adWt, Xq = xq, 



(20) 




Y = x - a, 



(21) 



ne Kt (X t 
ae Kt dW t 



a)dt + e Kt dX t 



It follows that {Yt,t > 0} is a continuous Gaussian process with variance function 





Wang and Potzelberger: Boundary Crossing Probabilities for Diffusion Processes 



9 



For every s > 0, define W s := Y t r s \ — Yq. Then, it follows from Wo = and 

V(W s ) = V(Y t(s) ) = S (t( S )) = s, 

that {Ws, s > 0} is a standard BM. Furthermore, for every t > 0, Y t = Yq + W 8 u) and by (f2Tj) we 
have 

X t = e- Kt Y t + a 

= e' Kt (x -a + W s(4) ) +a. 

Therefore, we can write the one-sided BCP for {Xt} as 

P(X t < b(t),Vt G [0,T]) = P[e- Kt (x Q -a + W s{t) ) +a<6(t),Vt€ [0,T]) 

= P (W s(t) < a - x + e Kt (6(t) - a), Vt € [0, T]) 
= P (W s < a - x + e"*W (&(*(*)) - a), Vs G [0, S]) , 

where 5 = s(T) = a 2 (e 2KT - l) /2 k. Obviously, similar relation for the two-sided BCP can be 
established in the same way. Thus, we have the following result. 

Corollary 7 Let {X t ,t > 0} be an O-U process satisfying (2(h) . Then there exists a standard BM 
{W s ,s > 0}, such that for any boundaries a(t),b(t), 

P x (a,b,T)=P w (c,d,S), (22) 

where 

c(s) = a — xq + [a(t(s)) — a] ll + j , (23) 

/ 2ks\ 1/2 

d(s) = a-x + [b(t(s)) - a ]ll + — \ , (24) 

t( a ) = J-log(l + ^V a >0 (25) 



2k"°\ ' a 2 , 
and S = a 2 (e 2nT - l) /2k. 

The above derivation can be extended to more general O-U processes with time-dependent 
coefficients, i.e., processes which satisfy 

dX t = K(t)(a(t) - X t )dt + a(t)dW t , X = x , (26) 

where n(t),a(t) : [0, oo) i— > iR+ and a(t) : [0, oo) i— ► iR are real, deterministic functions. Indeed, if 
we define, for every i > 0, 

Y t = exp Qf (X t - j(t)), Y = x - 7(0), 

where 

7 (t) = e ~'a(0) + e~* / e u a(u)du, 
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then by Ito's lemma, we have 

dY t = exp (^J K(v)dvJ a(t)dW t . 

It follows that 

s(t) = V(Y t ) = £exp (2 J\{v)dv^j a 2 (u)du, (27) 
which is strictly increasing and hence admits the inverse function t(s),s > 0. Therefore, 
X t = 7(t) + exp (- J K (u)du \ Y t 

= + exp (- J K(u)d\ (x - 7(0) + W s(t) ) , 

where W s = Y t / S ) — Yq is a standard BM. Thus, we have the following result. 

Corollary 8 Let {Xt,t > 0} be an O-U process satisfying )126\). Then there exists a standard BM 
{W s ,s > 0}, such that for any boundaries a(t),b(t), Px(a,b,T) = Py^(c,d, S), where 

( f t{s) \ 
c(s) = 7(0) - x + [a(t(s)) - j(t(s))] exp / n(u)du , 



d(s) = 7(0) - x + [&(*(«)) - 7(i(a))] exp I /* ( * As(t*)d« j , 



is t/ie inverse function of s(t) defined in (21) and S = s(T). 

Relation (|22|) allows us to derive exact formulas for the BCP of an O-U process {Xt} crossing 
ceratin boundaries, using the existing results for Brownian motion crossing, e.g., linear or square- 
root boundaries. For the simplicity of notation, in the following we present three examples for 
one-sided BCP. 

Example 9 First, consider boundaries a(t) = —00 and b(t) = a + he Kt ,0 <t<T, where h £ M is 
an arbitrary constant. Then by (I23h and (I24p . c(s) = —00 and d(s) = a — XQ + h+2hns/a 2 is a linear 
function in s G [0, S], where S = o~ 2 (e 2KT — 1)/2k. Therefore, by (|22p and the well-known formula 
for the BCP for BM crossing a linear boundary (e.g., equation (3) of Wang and Potzelberger 1997), 
we have 

p(hT , ' he 2 ^ + a-x \ ( 4hn(h + a- X0 ^ Jhe^-a + X0 -2h 
P x {a, b,T) = $ I ? == = I - exp ( ) $ 



a^(e 2KT - 1)/2k) V ^ / y ^/(e 2 ^ - 1)/2k 

where $ is the standard normal distribution function. 

Example 10 Another set of boundaries are a{t) = —00 and b(t) = a + he~ Kt ,0 <t<T. In this 
case c(s) = —00 and d(s) = a — xq + h is a constant. Therefore, again by (|22p we have 



Px{a,b,T) =2$ 



a — xq + /i 



<7J(e 2KT -l)/2 K/ 
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Example 11 Further, consider constant boundaries a(t) = — oo and b(t) = h, < t < T. Then 

c(s) = — oo and 

/ 2ks\ 1 I 2 
d{s) = a- xo + (h — a) ( 1 + \ 

is a square-root boundary, for which the "exact" BCP for BM is known (Daniels 1996). A numerical 
example of this case is given in section 7. 



5 Growth processes 

Another important stochastic model in population genetics is the growth processes (Ricciardi et al 
1999), which is defined by the SDE 

dX t = (aX t - (3X t log X t )dt + aX t dW t , X = x (28) 

and has state space H+, where a, (3 and a are positive constants. It is easy to see that the growth 
processes belong to the G-class described in Corollary [5j 

Similar to the O-U process, again by equations (fT5|) and (fT6|) . functions f(t,x) and at can be 
determined as 

f(t,x) = — [logx + ^-j 

and at = e"*. Hence the time transformation is 

s (t) = f e 2/3u du = -*^(e 2/3t - 1), t > 0. 
Jo 2(3 

Therefore, by (jlip and (|12p . we have the following result. 



Corollary 12 Let {Xt,t > 0} be a growth process satisfying \28\) . Then there exists a standard 
BM {W s ,s > 0}, such that for any boundaries a(t),b(t), Px(a,b,T) = P^(c,d, S), where 

c(s ) = V1 + 2/js ^log a(t(s)) + ) - ~ [ log x + "-^f^ ) , (29) 

d{s) = V1 + 2/3S L gb(t(fl)) + I -I log,,, + - ) . C!0) 



(7 



o- 2 - 2a\ 


\~\\ 


^log 


+ 


a 2 - 2a 


2(3 j 




2/3 


a 2 -2a\ 


\-\\ 


log 


S^o + 


cr 2 -2a 


2/3 j 




2(3 



*(s) = ^log(l + 2/fr) (31) 



and 5 = (e 2 ^ T - l)/2/3. 

Now we derive explicit formulas for the BCP for some special boundaries. 
Example 13 Let a(t) = and 

6(f) = «p L^-^^Vo<i<T, 

where /j is a constant. Then by (|29|) and ([30]) . c(s) = — oo and 

,, N 2h(3s if,, a 2 - 2a\ n ^ ^ _ 
d(s) = —£- + - U-logxo ^— j ,0 < s < S 
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which is a linear function and S = {e 2 ^ T — l)/2(3. Therefore, 
'2(3(he^ T -logx )-a 2 + 2a 

P x {a,b,T) = 3> 



a^2[3(eW T - 1) 

^ I ih£Qog x -h) + 2h(a 2 - 2a) \ $ ^ 23(he 2l3T -2h + log gg) + q- 2 - 2a 



a' 



\ aJ2(3{e 2 P T - 1) 



where $ is the standard normal distribution function. 
Example 14 Let a(t) = and 

b(t) = exp (he'^ - - ~ | , < t < T, 



213 

where h is a constant. Then c(s) = — oo and 



d(s) = \{h- log - a \p a ) , < s < 5. 



ft(a,6,T) = 2d> I 2 ^- ; lQg ^-" 2 + 2a I - 1. 



Therefore, 

Px(a,&,T) = 2<D| M 

a^2fi{e 2 ^ T - 1) 

Example 15 Again, consider the constant boundary b{t) = h > 0,0 < t < T. Then 

= [ l0gh + ^23—) " a { lOgX0 + ^T) '° " S " S 

with S = (e 2/3T — l)/2/3. A numerical example of this case is given in section 7. 

6 Geometric Brownian motion 

One popular stochastic model in mathematical finance is a generalization of the classical Black- 
Scholes model (Black and Scholes 1973) to the time-dependent interest rate process, as defined 
by 

dX t = r{t)X t dt + aX t dW t , X = x , (32) 

where r(t) : [0, oo) i— > iR+. It is well-known that in this case there exists a risk-neutral probability 
measure, under which the underlying process is a geometric BM 

X t = X exp \R(t) - ^ + <?w)j , 

where R(t) = ^r{u)du and {Wt} is a standard BM under the same risk-neutral probability mea- 
sure. In the following we use the approach of this paper to derive the BCP for {A^}. 
First, since equation (|16p now becomes 

9/ = at_ 
dx ax ' 
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function f(t,x) is given by 



f(t,x) = ^ (log x + h(t)) 
a 



where h(t) is an arbitrary function to be determined later. Further, since 

df = (log x + h(t)) da t a t dh(t) 
dt a dt a dt ' 

and 



d 2 f 



dx 2 ax 2 ' 

equation (fT5|) becomes 



which holds, if at = 1 and 

dh(t) _ a^ 
dt ~ T " 

The last equation above has a particular solution 

t f a 2 



r(t). 



h(t) = ^ y--r(u)jdu 

- q - m 

Therefore, we obtain at = 1 and 

f( t , X )=i(iog X +^-R(t)y 

By Theorem[2]the time transformation is s(t) = t, t > 0, and by (jlip and (|12p . we have the following 
well-known result. 

Corollary 16 Let {Xt,t > 0} 6e a geometric BM satisfying A32\). Then there exists a standard 
BM {W s , s > 0}, such that for any boundaries a(t),b(t), Px(a,b,T) = Py^(c,d,T), where 

c(t) = 1 (log (^) + ^ - , < t < T, (33) 

and 

d(t) = - I log ( ^ ) + ^— - fl(t) I ,0 < i < T. (34) 



^ og (m )+ q. m y 0< _ t< _ T . 



Now we consider some special cases. 

Example 17 Let a(t) = and b(t) = exp(pt + q + R(t)), where p, q are constants. Then by ([33 
and PI]) , c(t) = — oo and 

= - fp+ + ^i^io,o < t < T. 

a \ 2 1 a 
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Therefore, 

' (p + o 2 /2)T + q-\ogx Q 



P X (a,b,T) = $ 



( (2p + a 2 ) (log Xo - q ) \f (p + a */2)T - q + log x \ 
exp l * 2 ) { o-VT J' 

Example 18 Another special case is the constant interest rate r(t) = r, so that R{t) = rt. In this 
case, if we take constant boundaries a(t) = and b(t) = h > 0, then c(t) = — oo and 

d W = K^" r ) t+ - l0S (^)'°- t - T ' 

Therefore, 

\a 2 /2-r)T + \og{h/x )\ 



P X (a,b,T) = $ 



(2r-a 2 )log(h/x ) \ ( (a 2 /2 - r)T - \og(h/x Q ) \ 



exp . 



More numerical examples for the geometric and standard BM are given in the next section. 



7 Numerical examples 

In this section we compute some numerical examples for the BCP for various processes discussed 
in previous sections. For the simplicity of notation, we present examples of one-sided BCP only. 

Given any diffusion process {Xt} and a boundary, the BCP for {Xt} can be given by the BCP 
for Brownian motion and a transformed boundary b, i.e., Pw(— oo, b, T). As discussed in section 
2, the BCP P\v(— oo, b, T) can be approximated by Pw{— oo, b n , T), where b n is a piecewise linear 
boundary converging to b uniformly on the interval [0, T]. In particular, if b n (t) and b n (t) are such 
piecewise linear functions approaching b(t) from below and above respectively, then 

P w (-oo,b n ,T) < P w (-oo,b,T) < P w (-oo,b n ,T). 

The two sides of the above inequalities are given by ([2]) with function g given in ([3|). The expectation 
in (j2J) can be easily computed through Monte Carlo simulation, because Wt x , Wt 2 , Wt n have a 
multivariate normal distribution. 

For all examples in this section, b n (t) and b n (t) are constructed with n = 128 equally-spaced 
nodes and lower and upper bounds for the corresponding BCP are computed. In each Monte Carlo 
simulation, 10 6 repetitions are carried out. 

Example 19 (O-U process) Consider constant boundaries of Example [TT| where a(t) = —oo 
and b(t) = h, < t < T. Then, c(s) = — oo and d(s) is a square-root boundary. In particular, 
we set parameters to a = xq, a 2 = 2k = 1, h = xq + 1 and T = 1. Then d(s) = \f\ + s, for 
0<s<5 = e — 1. Using the procedure described above, the computed lower and upper bounds 
for the BCP are 0.721463 < P x (a,b,T) < 0.721464, with a simulation standard error of 0.000440. 
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Example 20 (Growth process) Now consider boundaries of Example [T5l where b{t) = h, < 
t < T. Then d(s) is again a square-root boundary. In particular, if we take a 2 = 2a, (3 = 1/2, 
xo = 1, h = e a and T = 1, then d(s) = VI + s and 5 = 6 — 1, which is the same as the square-root 
boundary of Example 20. The lower and upper bounds for the BCP are therefore given there. 

Example 21 (Geometric BM) Roberts and Shortland (1997) give an example of pricing a Eu- 
ropean call option, who's underlying security Xt follows the SDE (|32p with a = 0.1, r(t) = 
0.1 + 0.05e _t and Xq = 10. The option has a knock-in boundary b(t) = 12 and maturity T = 1. 
Then by Corollary [161 and equation (|34"|) . we have Px(—oo,b,T) = P^(- oo, d, T), where 



Therefore, using the approach of this paper the lower and upper bounds for the BCP are found to 
be 0.603728 < P x (-oo,b,T) < 0.603729, with a simulation standard error 0.000483. 

Example 22 (BM and Daniels boundary) Finally, consider the well-known Daniels' bound- 
ary for standard BM W = {W t ,t > 0} (Daniels 1969, 1996), which is 



with T = 1. This boundary has been used by many authors to test their computational algorithms. 
The "exact" BCP for this boundary is known to be Pw {— oo, b, T) = 0.520251. The procedure 
described above gives an approximation for the corresponding BCP as Pw(—oo,b,T) 0.520293 
with a simulation standard error 0.000490. 

8 Conclusions and Discussion 

We proposed a direct method for computing boundary crossing probabilities (BCP) for a class of 
diffusion processes which can be expressed as piecewise monotone functionals of a standard Brow- 
nian motion. This class includes many interesting processes in real applications, e.g., Ornstein- 
Uhlenbeck, growth processes and geometric Brownian motion with time dependent drift. This 
method applies to both one-sided and two-sided general nonlinear boundaries, which may be dis- 
continuous. Using this approach explicit formulas for boundary crossing probability for certain 
nonlinear boundaries have been obtained, which are useful to evaluate and compare computational 
algorithms. 

Compared to the traditional approach for the first-passage time densities, the BCP method can 
be applied to both one-sided and two-sided boundaries. It also applies to discontinuous bound- 
aries. The actual computation can be easily done using Monte Carlo simulation methods, which 
is practical and easy to implement. Furthermore, the approximation errors are also computed 
automatically. The approach of this paper is based on the stochastic differential equation for the 
processes. It is different from the traditional methods of first-passage time densities, which is based 
on the partial differential equation for the transition density function. 
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d(t) = 101og(1.2) - 0.5 - 0.95£ + 0.5e~',0 <t<T. 
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